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Density of states, dynamic (optical) conductivity and phase diagram of paramagnetic two - di- 
mensional Anderson - Hubbard model with strong correlations and disorder are analyzed within 
the generalized dynamical mean field theory (DMFT+S approximation). Strong correlations are 
accounted by DMFT, while disorder is taken into account via the appropriate generalization of the 
self - consistent theory of localization. We consider the two - dimensional system with the rect- 
angular "bare" density of states (DOS). The DMFT effective single impurity problem is solved by 
numerical renormalization group (NRG). Phases of "correlated metal", Mott insulator and corre- 
lated Anderson insulator are identified from the evolution of density of states, optical conductivity 
and localization length, demonstrating both Mott - Hubbard and Anderson metal - insulator tran- 
sitions in two - dimensional systems of the finite size, allowing us to construct the complete zero 
- temperature phase diagram of paramagnetic Anderson - Hubbard model. Localization length in 
our approximation is practically independent of the strength of Hubbard correlations. However, the 
divergence of localization length in finite size two - dimensional system at small disorder signifies 
the existence of an effective Anderson transition. 

PACS numbers: 71.10.Fd, 71.27. +a, 71.30.+h 



I. INTRODUCTION 

The study of disordered electronic systems with the 
account of interaction effects belongs to the central prob- 
lems of the modern condensed matter theory^. According 
to the scaling theory of localization- there is no metal- 
lic state in two - dimensional (2D) systems, all electronic 
states are localized at the smallest possible disorder. This 
prediction was made first for noninteracting 2D systems, 
soon after it was shown that the weak electron - electron 
interaction in most cases enhance localization 3 . Experi- 
ments performed in early 80 's on different 2D systems 4 
essentially confirmed these predictions. However, some 
theoretical works produced an evidence of rather differ- 
ent possibility 5 for 2D systems to evolve to the state 
with even infinite metallic - like conductivity at zero 
temperature (T = 0) in case of weak disorder and suffi- 
ciently strong correlations. Experimental observation of 
metal - insulator transition (MIT) in 2D systems with 
weak enough disorder but strong correlations (low elec- 
tronic densities) 6 , which apparently contradicted the pre- 
dictions of the scaling theory of localization, stimulated 
extensive theoretical studies with no widely accepted so- 
lution up to now (see the review in Ref. 0). 

One of the basic models allowing for the simultane- 
ous account of both strong enough electronic correla- 
tions, leading to Mott MIT transition 8 -, and effects of 
strong disorder, leading to Anderson MIT—, is the An- 
derson - Hubbard model, intensively studied in recent 



In Refs! 10 i n ' 12 three - dimensional (3D) Anderson - 
Hubbard model was analyzed with dynamical mean field 
theory (DMFT), which is extensively used in the theory 
of strongly correlated electron o 17 i 18 ' 19 i 20 . However, disor- 
der effects were mostly taken into account via the average 
density of states and the coherent potential approxima- 
tion (CPA) 21 i 22 , which misses the effects of Anderson lo- 



calization. To overcome this problem Dobrosavljevic and 
Kotliar— have proposed the DMFT approach, where the 
solution of self - consistent stochastic DMFT equations 
were used to calculate geometrically averaged local den- 
sity of states. This approach was further developed in 
Refsj U i 12 with DMFT account for Hubbard correlations, 
which leads to highly nontrivial phase diagram of 3D 
paramagnetic Anderson - Hubbard model 12 , containing 
the phases of correlated metal, Mott insulator and corre- 
lated Anderson insulator. However, the major problem 
of the approach of Refsj 10 ' 11 ! 12 is its practical inability of 
direct calculations of conductivity, which actually deter- 
mines MIT itself. 

In the previous worki^ we have studied the 3D param- 
agnetic Anderson - Hubbard model using our recently de- 
veloped DMFT+S approximatio n 23 ' 24 ^^, which con- 
serves the standard single impurity DMFT approach, 
taking into account the local Hubbard correlations, allow- 
ing the use the usual " impurity solvers" like NRG — i 28 ' 29 , 
at the same time allowing to include additional (local or 
nonlocal) interactions. Strong disorder was accounted 
via some generalization of the self - consistent theory of 
localizatio n 30 ' 31 ! 32 ' 33 ! 34 ' 35 . In the framework of this ap- 
proach we have been able not only to reproduce the phase 
diagram qualitatively similar to that obtained in Ref. IT3 , 
but also calculate the dynamic (optical) conductivity for 
the wide frequency range. 

In Ref. Il5| the Hubbard - Anderson model was stud- 
ied both for 3D and 2D systems. As the main mecha- 
nism leading to derealization a kind of "screening" of 
the random (disorder) potential by local Hubbard in- 
teraction was introduced^. Then the Anderson - Hub- 
bard model was reduced to an effective single - particle 
Anderson model with renormalized distribution of local 
site energies, which was calculated in the atomic limit. 
All the other effects of electronic correlations were ne- 
glected. Strong disorder effects were accounted within 
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self - consistent theory of localization. In this approach 
the authors obtained the significant growth of localiza- 
tion length with growing Hubbard interaction in 2D. 
However, localization length itself remained finite, the 
system being localized at smallest possible disorder, so 
that Anderson in 2D is still absent. Similar result was 
obtained also in numerical simulations of 2D Anderson - 
Hubbard model in Ref. Il6l . 

In this work we present a direct generalization of the 
method used in Ref. [HI to the case of 2D systems. We 
use the DMFT+E approach to calculate DOS, optical 
conductivity, localization length and construct the phase 
diagram of 2D paramagnetic Anderson - Hubbard model 
with strong electronic correlations and strong disorder. 
Strong correlations are taken into account via DMFT, 
while disorder effects are treated by the appropriate gen- 
eralization of the self - consistent theory localization. 

The paper is organized as follows: in section [H] we 
present a brief description of our DMFT+E approxima- 
tion as applied to disordered Hubbard model. In section 
IIIII we formulate the basic DMFT+S expressions for op- 
tical conductivity and self - consistency equation for the 
generalizes diffusion coefficient. Our results for DOS, op- 
tical conductivity and localization length are given in sec- 
tion IIV[ where we also analyze the phase diagram of 2D 
disordered Hubbard model and briefly discuss the optical 
sum rule within our approach. Finally we present a short 
conclusion, which includes the discussion of problems yet 
to be solved. 



II. BASICS OF DMFT+E APPROACH. 

In the following we consider paramagnetic disordered 
Anderson - Hubbard model at half - filling for arbitrary 
correlations and disorder. This model treats both Mott 
- Hubbard and Anderson MIT on the same footing. The 
Hamiltonian of the model can be written as: 

H = -t a\ a a ia + ^2 e i n i° + ^ E n ^ n ii> W 

(ij)<7 i<y i 

where t > is nearest neighbor transfer integral, U is the 
local Hubbard repulsion, rii a = a^aia is electron num- 
ber operator at a given site i, (aJ CT ) is annihilation 
(creation) operator for an electron with spin a, local en- 
ergies £j are assumed to be randomly and independently 
distributed on different lattice sites. To simplify diagram 
technique in the following we assume distribution to 
be Gaussian: 

Here A serves as disorder parameter and Gaussian ran- 
dom field ("white" noise) of energy levels q at different 
lattice sites induces "impurity" - like scattering, leading 
to the standard diagram technique for calculations of the 
averaged Green's functions^ 5 .. 



DMDF+E approach, initially proposed in 
Refsj 23 i 24 ' 25 i 26 as a simple method to include non- 
local interactions (fluctuations) into the standard (local) 
DMFT scheme, is also very convenient for the account in 
DMFT of any additional interaction (local or nonlocal) 
of arbitrary nature. 

In DMFT+S approximation we choose the single - 
particle Green's function in the following form: 

Gp(£) = e + /i-e(p)-S(e)-Sp( £ )' (3) 

where e(p) is the "bare" electron spectrum, E(e) is the 
local (DMFT) self - energy due to Hubbard interactions, 
while E p (e) is an "external" (in general case momentum 
dependent) self - energy due to some other interaction. 
The main assumption of our approach (both its advan- 
tage and deficiency) is precisely in this additive form (ne- 
glect of interference effects) for the total self - energy in 
(j3j23,24,25.26^ which allows us to conserve the standard 
form of self - consistent DMFT equations 2 ^ with two 
major generalizations. First of all, at each iteration of 
DMFT - loop we recalculate an "external" self - energy 
E p (/i, £, [G p (e)]) within some (approximate) scheme, tak- 
ing into account the "external" interaction (in the present 
case that due to disorder scattering). Secondly, the lo- 
cal Green's function for an effective DMFT - impurity 
problem is defined as: 

G «( e ) = ^E e+M _ e(p) _ s(e) _ Ep(e) ' ( 4 ) 

at each step of the standard DMFT procedure. Finally 
we obtain the desired Green's function in the form of Eq. 
([3]), where E(e) and E p (e) are self - energies obtained at 
the end of our iteration procedure. 

For S p (e) appearing due to disorder scattering we shall 
use the simple one - loop contribution, neglecting dia- 
grams with "crossing" interaction lines, i.e. self - consis- 
tent Born approximation 35 , which in the case of Gaussian 
disorder ^ leads to the usual expression: 

E p (e) = A 2 G{e, p) = E imp (e) (5) 
p 

and "external" self - energy in this case is p-independent 
(local). 



III. OPTICAL CONDUCTIVITY IN DMFT+E 
APPROACH. 

It is obvious that calculations of optical (dynamic) con- 
ductivity provide the direct way to study MIT as fre- 
quency dependence of conductivity, as well as its static 
value at zero frequency of an external field, allows the 
clear distinction between metallic and insulating phases 
(at temperature T = ). 
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FIG. 1: Diagrammatic representation of $® RA (uj,q) 



Local nature of irreducible self - energy in DMFT al- 
lows to reduce the problem of calculation of optical con- 
ductivity to calculation of the usual particle - hole loop 
without DMFT vertex corrections due to local Hubbard 
interactio n 13 ! 26 . The final expression for the real part of 
optical conductivity obtained in this way in Refsi 13 ' 26 , 
takes the following form: 



Rea(u)) 
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e u> 
~2T 



xRe 4> ora (uj) 



de [/(£_) -/(£+)] x 



E*(e + ) - E«(e_) 



(6) 



Here f(e) is Fermi distribution, e± = e ± £ and 



^O-Rfi(itVL), \ A 0i?,ft(i?A), 

$e ; (w,q)-$ e (W,0) 



where we have introduced two - particle loops: 



(7) 



p 

r^(^)( £ _,p_; £+) p + ), (8) 

represented diagrammatically in Fig. [1] with p± = 
p ± 3 and i? and A indices corresponding to re- 
tarded and advanced Green's functions. The vertices 
r- RR (' RA )(£_ ) p_;£ + ,p_|_) contain all vertex corrections 
due to disorder scattering, but do not include vertex cor- 
rections due to Hubbard interaction. 

Thus, the problem is much simplified. To calculate op- 
tical conductivity in DMFT+E approximation we have 
only to solve single - particle problem to determine the 
local self - energy E(£±) with the help of DMFT+E pro- 
cedure, described above, while the nontrivial contribu- 
tion of disorder scattering enters via (fp RR ( RA ) of Eq. ([7]) , 
which may be calculated in some appropriate approxima- 
tion. In fact, (j) 0RR (. RA ) contain only disorder scattering, 



though using as the "bare" the Green's functions includ- 
ing the DMFT self - energies, already determined with 
the help of DMFT+E procedure. Eq. ([6|) guarantees the 
effective interpolation between the case of strong corre- 
lations in the absence of disorder and the case of pure 
disorder in the absence of Hubbard correlations. 

The most important &® RA (uj 7 q) loop may be calcu- 
lated using the basic approach of the self - consistent 
theory of localizatio n 30 i 31 ' 32 i 33 i 34 i 35 with some general- 
izations accounting for Hubbard interaction within the 
DMFT+E approach^ 3 .. 

The rest is the direct generalization of the scheme pro- 
posed in Ref. [l3| for the two - dimensional case. Here we 
present only some basic points of the approach of Ref. [H, 
stressing important differences due to two - dimension- 
ality of the model. 

In RA-channel the two - particle loop Q® RA (q,(li) pos- 
sesses a diffusion - like contribution: 



<J>r A (q,£) 



uj + iD(uj)q 2 ' 



(9) 



where AG P = G R (e+,p) - G A (e-,p). The important 
difference with noninteracting case is contained in 

£ = £+-£_-£*(£+) = 

cj-E fi (e + ) + E A (e_) = w- AE KA (w) (10) 

which replaces the usual w-term in the denominator of 
the standard expression for $>® RA (ui,q)2^. 
Then ^ can be rewritten as: 



xRe ■ 
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(11) 



The second term in (jllj) . which is actually irrelevant at 
small uj, can be obtained from ([7]) calculating $® rr (uj, q) 
in the usual "ladder" approximation. 

Repeating the derivation scheme of the self - consis- 
tent theory of localization presented in Ref. fl3l . we obtain 
the following equation for the generalized diffusion coef- 
ficient: 



D(u) 



+ A^(AG P ) 2 £, 



uj + iD(uj)q 2 



(12) 



where d = 2 is spatial dimension and AY, RA p (. 



Y> R mp (e + ) — E^ p (£_) is determined by disorder scatter- 
ing. The average velocity < v >, well approximated by 
Fermi velocity, is defined as: 



< v >- 



Ep|v P |AGp 
E p AGp ' 



de(p) 
dp 



(13) 
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Due to the limits of diffusion approximation summation 
over q in (Tl2"|) should be limited b y 33 i 35 : 



q < fco = Min{l 1 ,pf} 



(14) 



where I =< v > /2*y is the mean - free path due to elastic 
scattering (7 is the scattering rate due to disorder), pp is 
Fermi momentum. In our two - dimensional model An- 
derson localization takes place at infinetisimal disorder. 
However, for small disorder localization length is expo- 
nentially large, so that the size of the sample becomes 
important. The sample size L may be introduced into 
the self - consistent theory of localization as a cutoff of 



diffusion pole contribution at small q du i , i.e. for 
q~k L = 1/L. 



(15) 



Eq. (TT!?]) for the generalized diffusion coefficient re- 
duces to a transcendental equation, which is easily solved 
by iterations for each value of u, taking into account that 
for d = 2 and cutoffs defined by Eqs. (jT4"|) . (fl~5j) the sum 
over q in (fT2|) reduces to: 



E 

q 



1 



1 



u) + iD(uj)q 2 i2ttD(ijj) JhL_ y 2 + 



ydy 



(16) 
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Solving Eq. (fl2|) for different values of our model pa- 
rameters and using Eq. (jTTJ) we can directly calculate 
the optical (dynamic) conductivity in different phases of 
the Anderson - Hubbard model. 

For uj — ► (and at the Fermi level (e — 0) obvi- 
ously also Co — > 0), in Anderson insulator phase we ob- 
tain the localization behavior of the generalized diffusion 
coefficient^M^: 



D(ui) = -iCbRio 



(17) 



After the substitution of (TT7|) into Q12[) we get an equa- 
tion, determining the localization length Ri oc : 



Ric 



<v>^ 
dA 4 



. p 



(AG P ) 2 E r 



1 



2 a 2 



Rio 

(18) 

where the sum over q is defined by (|16p . As we shall see 
in the following, for an infinite two - dimensional system 
(L — > 00) localization length, determined by Eq. [)18p 
remains finite (though exponentially large) for the small- 
est possible disorder, signifying the absence of Anderson 
transition. However, for the finite size systems localiza- 
tion length diverges at some critical disorder, determined 
for each value of the system size L. Qualitatively, this 
critical disorder is determined from the condition of local- 
ization length (in infinite system) becoming of the order 
of characteristic sample size 

R fo7°° ~ L - Thus > in nnitc 
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FIG. 2: Density of states of the half - filled Anderson - Hub- 
bard model for different values of U and A = (no disorder). 



two - dimensional systems the Anderson transition and 
metallic phase do exist for disorder strength lower, than 
this critical disorder. In the following, this kind of metal- 
lic phase will be referred to as a phase of the "correlated 
metal" in finite 2D systems. 



IV. MAIN RESULTS 

Below we present the results of extensive numerical 
calculations for 2D Anderson - Hubbard model on the 
square lattice with rectangular "bare" DOS, correspond- 
ing to the bandwidth W = 2D: 



N (e) = 



i 

2D 




|e|<D 
lei > D 



(19) 



The choice of this model DOS is dictated by its 2D na- 
ture. 

Everywhere below we give the values of DOS in units 
of the number of states per energy interval, per lattice 
cell of the volume a 2 (a - lattice parameter), per one 
spin projection. 

As we concentrate on half - filled case, Fermi level is 
always assumed to be at zero energy. 

As "impurity solver" for DMFT effective impurity 
problem we have used the reliable numerical renormal- 
ization group (NRG) approac h 27 ! 28 ' 29 . Calculations were 
made for low enough temperature T ~ 0.001D, so that 
temperature effects in DOS and conductivity are just 
negligible. 

Now we present only most typical results. 
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FIG. 3: Density of states of the half - filled Anderson - Hub- 
bard model for different values of disorder A and U/2D = 
1.25, typical for correlated metal (in the absence of disorder). 



A. Evolution of the density of states 



0.8 
0.6 
0.4 
0.2 







1 A/2D=0 

2 0.1 1 

. 3 0.19 

4 0.25 

5 0.37 

6 0.5 




I'i U/2D=2 


M 













J-p-l — I — . — I — . — I — . — I > I f I — . — I — . — I — I — I — ^-1 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

e/2D 



FIG. 4: Density of states of the half - filled Anderson - Hub- 
bard model for different values of disorder A and U/2D — 2, 
typical for Mott insulator (in the absence of disorder). At 
the insert - restoration of quasiparticle band by disorder in 
coexistence (hysteresis) region for U = 1.5D, obtained from 
Mott insulator with decreasing U. 



In Fig. [5] we show evolution of the DOS with the 
growth of Hubbard interaction U in the absence of dis- 
order. At small U (curve 1 in Fig[2]) we observe practi- 
cally rectangular DOS almost coinciding with the "bare" 
one. As U grows a typical three peak structure of DOS 
appear o 19 ' 20 i 29 (curves 3,4,5 in FigJ5]): a narrow quasi- 
particle peak at the Fermi level with upper and lower 
Hubbard bands at s ~ ±U/2. Quasiparticle peak nar- 
rows as U grows in metallic phase, disappearing at Mott 
MIT for U = U c2 ps 1.83W. With further growth of U 
(curves 6,7 in Fig(2]) dielectric gap opens at the Fermi 
level. 

In Fig. [3] we show the results for DOS obtained at rel- 
atively weak correlation strength U = 1.25W (W = 2D), 
so that the system is rather far fron the Mott transi- 
tion, but for the wide range of disorder strength A. We 
observe typical widening of the band with appropriate 
suppression of DOS as disorder grows. 

In Fig. 2] DOS evolution is shown as disorder A grows 
at U = 2W, typical for Mott insulator in the absence of 
disorder. It can be seen that the growth of disorder leads 
to restoration of the quasiparticle peak in DOS. Similar 
unusual behavior of DOS (closing of dielectric gap by 
disorder) was first noted in 3D systems^. However, in 
the present 2D case it docs not, in general, signify the 
transition to the correlated metal phase, at least for the 
infinite systems we are, in fact, dealing with correlated 
Anderson insulator (cf. below). 

The physical reason for this unusual restoration of the 
quasiparticle peak in DOS is clear. Controlling parame- 
ter for appearance or disappearance of quasiparticle peak 
in DMFT in the absence of disorder is the ratio of Hub- 
bard interaction U and the "bare" bandwidth W = 2D. 



Disordering leads to the growth of the effective band- 
width W e f f (in the absence of Hubbard interaction) and 
appropriate suppression of U/W e ff ratio, which obvi- 
ously leads to the restoration of quasiparticle band in 
our model. In more details this qualitative picture will 
be discussed in Section IIVC1 where our results for DOS 
will be used in construction of the phase diagram of 2D 
Anderson - Hubbard model. 

It is well known, that hysteresis behavior of DOS 
is obtained for Mott-Hubbard transition if we perform 
DMFT calculations with U decreasing from insulating 
phas o 20 i 29 . Mott insulator phase survives for the val- 
ues of U well inside the correlated metal phase, ob- 
tained with the increase of U. Metallic phase is re- 
stored at U c i ps 1A2W. The values of U from the in- 
terval U c i < U < U C 2 are usually considered as belong- 
ing to coexistence region of metallic and (Mott) insulat- 
ing phases, with metallic phase being thermodynamically 
more stabl o 20 i 29 . In the coexistence region disorder in- 
crease also leads to the restoration of quasiparticle peak 
in the DOS (see insert of Figf4]). 



B. Optical conductivity: Mott — Hubbard and 
Anderson transitions. 

The real part of optical conductivity was calculated 
for different combinations of parameters of the model, 
directly from Eqs. (fTTj) and (fT2"]) . using the results of 
DMFT+S procedure for single - particle characteristics. 
The values of conductivity below are given in natural 
units of e 2 /H. 

In the absence of disorder we just reproduce the results 
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FIG. 5: Real part of dynamic conductivity of the half - filled 
Anderson - Hubbard model for different values of disorder A 
and U/2D = 1.25, typical for correlated metal (in the ab- 
sence of disorder). At the insert - the same, but in wider 
frequency range. Thin dashed lines represent the results of 
ladder approximation. 



of the standard DMFT with optical conductivity charac- 
terized by the usual Drude peak at zero frequency and a 
wide maximum at uj ~ U, corresponding to transitions 
to the upper Hubbard band. As U grows Drude peak is 
suppressed and disappear at Mott MIT, when only re- 
maining contribution is due to transitions through the 
Mott - Hubbard gap. 

Introduction of disorder leads to qualitative change of 
frequency behavior of conductivity. Below we mainly 
present results obtained for the same values of U and 
A, which were used above to illustrate evolution of DOS. 

In Fig. [5]we show the real part of optical conductivity 
in 2D half - filled Anderson - Hubbard model for different 
disorder strengths A and U — 1.25W, when the system 
is far from Mott MIT. Thin dashed lines in Fig. [5] (as 
well as in the following figures) we show the results of the 
"ladder" approximation. In 2D model under considera- 
tion conductivity at zero frequency is always zero, and 
in contrast to 3D case^, even for the weakest possible 
disorder the peak in optical conductivity is at finite fre- 
quency. In the "ladder" approximation, which does not 
contain localization corrections, conductivity at to = is 
finite. Optical transitions to the upper Hubbard band at 
lu ~ U are practically unobservable in these data, only at 
the insert in Fig. [5j where we show the data for the wide 
frequency range, we can observe some weak maximum on 
curves I and 2, corresponding to transitions to the upper 
Hubbard band. 

In Fig. [6] we present the real part of optical conduc- 
tivity for different disorder strengths A and U — 2W, 
typical for Mott insulator. It can be seen from Fig. [6] 
that for small disorder we are in Mott insulator phase 



FIG. 6: Real part of dynamic conductivity of the half - filled 
Anderson - Hubbard model for different values of disorder 
A and U /2D — 2, typical for Mott insulator (in the absence 
of disorder). Curves 1,2 correspond to Mott insulator, while 
curves 3-6 correspond to correlated Anderson insulator. At 
the insert - localization behavior of conductivity. Thin dashed 
lines - ladder approximation. 



(curves 1,2), and with the growth of disorder in the ab- 
sence of Anderson localization (cf. thin lines correspond- 
ing to "ladder" approximation) we would be entering the 
metallic phase. However, in our model localization takes 
place at infinetesimal disorder and we are actually en- 
tering Anderson insulator phase, with conductivity go- 
ing to zero at zero frequency. Data in the frequency 
range corresponding to localization behavior of conduc- 
tivity a(u>) ~ to 2 are shown at the insert in Fig. [6] for 
curve 5 and 6, corresponding to large enough disorder. 
At small disorders the frequency region with localization 
behavior of conductivity is exponentially small 1 (which 
is due to the exponential growth of localization length at 
small disorder, cf. FigEjJ) and is practically unobservable. 

Dependence of optical conductivity on U is illustrated 
in Fig. [7j The growth of U shifts localization peak in con- 
ductivity to lower frequencies and leads to its narrowing. 
Apparently this is related to the appropriate suppression 
of quasiparticle peak width in DOS. The value of conduc- 
tivity at the maximum is U independent. It is interest- 
ing to note that for frequencies larger than the maximum 
position the growth of U suppresses conductivity, while 
for the frequencies lower than the maximum position the 
growth of U enhances conductivity (playing in a sense 
against localization). 

To confirm self - consistency of our approach to con- 
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TABLE I: Check of the single band optical sum rule in An- 
derson - Hubbard model. Optical integral is given in units of 

h 



A/2D 


7re 2 v~> & 2 £p 

— Z^p -a^- n p 


Wopt = J °° Rea(u>)cLj 


0.19 


0.099 


0.098 


0.25 


0.099 


0.098 


0.37 


0.092 


0.091 


0.5 


0.081 


0.082 



ductivity calculations we conclude this section with dis- 
cussion of the optical sum rule, which relates single - 
particle and two - particle characteristics^. 

The single - band Kubo sum rule2£ for dynamic con- 
ductivity can be written as: 

W op t = / Rea(Lu)dco = — V -jrfrip, (20) 

JO 1 p a Px 

where n p is single particle momentum distribution func- 
tion, determined by interacting retarded electron Green's 
function G R (e, p): 

1 f°° 

n p = / def(e)ImG R (e,p) (21) 

where /(e) is Fermi distribution. 

In Table I we show calculated values of the r.h.s. and of 
the l.h.s of Eq. (20]) for U = 1.5W. It is clearly seen that 
the optical (|20[) sum rule is fulfilled within our numerical 
accuracy. 

Very often the optical sum rule is understood as the 
equality of the optical integral W opt to the "universal" 

value of where uj p i is plasma frequency, which is 



FIG. 8: Dependence of the normalized optical integral of An- 
derson - Hubbard model on U for different values of disorder 
A. At the insert - similar dependence on A for different val- 
ues of U. Curves 1,2 - "correlated metal", transforming into 
Anderson insulator. Curves 3,4 - Mott insulator, obtained 
with the growth of U from "correlated metal" or Anderson 
insulator. 



strictly speaking is not correct in single band case and in 
this sense we may speak of optical sum rule "violation" . 
In fact the optical integral depends on parameters of the 
model, e.g. on Hubbard interaction and disorder (Fig. 
[8]| . The growth of U significantly suppresses the value of 
optical integral. Dependence on disorder strength A is 
also important, in particular disorder induced transition 
from Mott to Anderson insulator we observe a kind of 
discontinuity of optical integral (curves 3,4 at the insert 
in Fig. HJ. 



C. Localization length and phase diagram of 2D 
Anderson — Hubbard model. 

To proceed further on the left axis of Fig. [9]we present 
our data for the real part of conductivity at fixed and 
sufficiently small frequency ui — 0.00005D plotted as a 
function of disorder strength A. Circles show results of 
"ladder" approximation, triangles — results of self - con- 
sistent theory of localization (we take {7 = here) . Curve 
3, which practically coincides with results of the "ladder" 
approximation, is obtained from the usual Drude expres- 
sion: 

^)^(0)t^, (22) 
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FIG. 9: Left scale — dependence of conductivity on disorder 
A at fixed frequency u = 0.00005D and U/2D = 1. Circles 
- ladder approximation, triangles - self-consistent theory of 
localization. Curve 3, practically coinciding with the results 
of ladder approximation, is obtained from Drude expression 
(|22|) . Static conductivity for finite samples of sizes L = 10 8 a 
(curve 4) and L = 10 a (curve 5). Right scale — dependence 
of the logarithm of localization length on disorder A: for infi- 
nite sample (curve 1) and for finite samples of sizes L = 10 8 o 
(curve 2) and L = 10 5 a (curve 3). 

where the static conductivity is given by c(0) = 

e 2 N(0)D Q sa 4^T> with N (°) bein § the density of states 
at the Fermi level, Do is the classical diffusion coefficient, 
ep si D is Fermi energy. Impurity scattering rate was 
taken as 7 = 7rA^(0)A 2 sa ^jA 2 . It can be seen that 
the noticeable contribution of localization corrections to 
conductivity (clear difference between curve 2 curves 1 
and 3) appears only after conductivity drops below the 
values of the order of "minimal" metallic conductivity 
cr = 4- (our data for conductivity are actually normal- 
ized by this value in all figures) . We shall see below that 
precisely in this region a kind of Anderson MIT (diver- 
gence of localization length) takes place in 2D systems of 
reasonable finite sizes. 

On the right axis in Fig. [5] we show our data for the 
logarithm of localization length calculated from Eq. (fl~8"|) 
as a function of disorder for infinite sample (curve 1) 
and for finite samples with L — 10 8 a and L = 10 5 a 
(curves 2 and 3). It is clearly seen that localization length 
grows exponentially as disorder drops but remains finite 
in infinite 2D sample, signifying the absence of Ander- 
son transition. In finite samples localization length di- 
verges at some critical disorder (depending on the system 
size) demonstrating the existence of an effective Ander- 
son transition. From Fig. [5] it is seen that this critical 
disorder is achieved when localization length of an infi- 
nite system becomes comparable with characteristic size 
of the sample: Rf^ 00 ~ L. It should be noted that in 
our approach, opposite to the results of Ref. [IH, localiza- 



FIG. 10: Phase diagram of two - dimensional paramagnetic 
Anderson - Hubbard model at zero temperature. The bor- 
der of Mott insulator region f7 c 2(A) and the border of co- 
existence (hysteresis) region C/ c i(A) are obtained from the 
evolution of the density of states, [7* 2 (A) is calculated from 
(|24[) . Dashed stripe corresponds to the narrow region, where 
Anderson metal - insulator transition takes place in finite 
systems. 



tion length is practically independent of U, which leads 
to independence of critical disorder in 2D of correlation 
strength U . Similar result 2 was obtained in our approach 
for 3D systems^. 

On the left axis of Fig. [9] disorder dependence of static 
conductivity for finite samples of the sizes L = 10 8 a and 
L = 10 5 a (curves 4 and 5) is shown. For finite sys- 
tems with small disorder static conductivity is not zero 
(metal) . It gradualy goes down while disorder grows and 
becomes zero at the same critical value where localization 
radius diverges on the approach from insulating phase in 
a finite sample. Static conductivity of finite samples in 
our calculations practicaly does not depend on correla- 
tion strength U. Rather significant difference between 
the values of static conductivity and that of conductivity 
at small but non zero frquencies seen in Fig. [9] comes from 
exponential smallness of frequency range of localization 
behavior mentioned above. 

Let us now discuss our results for the phase diagram 
of 2D half - filled Anderson - Hubbard model, obtained 
from extensive DMFT+E calculations of DOS and anal- 
ysis of localization length behavior in finite 2D systems. 
The general form of this phase diagram in disorder - cor- 
relation (A, U) plane is shown in FigfTOl 

Dashed stripe in Fig. [TU] corresponds to the region of 



2 Calculations of localization length for 3D system performed by us 
after the publication of Ref. [13J have demonstrated its practical 
independence on the value of U 
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an effective transition from Anderson insulator to "metal- 
lic" phase. It boundaries were determined by divergence 
of localization length in finite samples with characteristic 
sizes L = I0 5 et (upper boundary) and L — f0 8 a (lower 
boundary) (cf. Figj9]). It should be stressed that further 
increase of system size, e.g. 10 times up to L = 10 9 a, 
leads only to practically negligible downshift of the lower 
boundary (decrease of critical disorder) of dashed stripe 
in Fig. \TU\ 

The dependence of i7 C 2(A), obtained from DOS behav- 
ior, determines the boundary for Mott transition and is 
defined by the disappearance of the quasiparticle peak in 
DOS and correlation gap opening at the Fermi level (cf. 
Figfi©. 

In our previous worki 3 - on 3D Anderson - Hub- 
bard model we have proposed a simple explanation of 
U c i tC 2(A) dependence. Assuming that the controlling pa- 
rameter of Mott - Hubbard transition given by the ratio 
of Hubbard interaction and effective bandwidth (depend- 
ing on disorder) ^'^^ is universal constant (indepen- 
dent of disorder), we get: 

Weff(A) W [ > 

where W e ff(A) is an effective bandwidth in the presence 
of disorder, calculated at U = in self - consistent Born 
approximation ([5]). In 3D modeU^ the dependence of crit- 
ical correlation strength on disorder i7 c i )C 2(A), obtained 
directly from the evolution of DOS, has shown quite satis- 
factory agreement with qualitative dependence obtained 
fromEq. (|2l5|) 3 

In 2D model under consideration here, solution of Eq. 
(25)) S ives: 

EGWA) = u c1 mo)^^ = 

= ^,c 2 (0)(^n(|±l) +C ), (24) 



where c = y 4 (A,) + 1. However, unlike in 3D casei 3 -, 
U C 2(A) dependence, obtained from DOS evolution is 
clearly different from the qualitative one U* 2 (A) (dotted 
line in Fig. [TO]) , determined by Eq. (|24)) . Probably, this 
is due to a significant change in the rectangular form of 
"bare" DOS with the growth of disorder A, which is ab- 
sent for semi - elliptic "bare" DOS used in 3D case in 
Ref. M 

As we already noted, with decrease of U from insu- 
lating phase Mott transition occurs at U — U c i(A) < 



3 Further extensive calculations performed after the completion of 
Ref.ll3l have confirmed practically ideal agreement between these 
dependencies. 



l7 c2 (A) and the coexistence (hysteresis) region is ob- 
served between U c i(A) and U C 2(A) curves on phase dia- 
gramm Fig lfOI 



V. CONCLUSION 

We have used the generalized DMFT+S approach to 
study basic properties of disordered and correlated An- 
derson - Hubbard model. Our method produces rela- 
tively simple interpolation scheme between two well stud- 
ied limits — that of strongly correlated Hubbard model 
in the absence of disorder (DMFT and Mott - Hubbard 
MIT) and the case of 2D Anderson insulator in the infi- 
nite system without electron - electron interactions. It 
seems that the proposed interpolation scheme reflects all 
the qualitative features of Anderson - Hubbard model, 
such as behavior of the density of states and dynamic 
conductivity. The general structure of the phase dia- 
gram obtained in DMFT+E approximation is also in rea- 
sonable agreement with the results of direct numerical 
simulations^. At the same time, DMFT+E approach is 
rather competitive in a sense of the amount of numer- 
ical work and allows direct calculations of all the basic 
observable characteristics of Anderson - Hubbard model. 

It should be stressed that an effective Anderson transi- 
tion obtained here in the case of finite size 2D systems is 
in no sense attributed to electronic correlations and fol- 
lows directly from self - consistent theory of localization 
also in the absence of correlations. 

The main shortcoming of the method used is the ne- 
glect of the interference between disorder scattering and 
Hubbard interaction, which leads to the independence 
of localization length and critical disorder A c (in finite 
2D systems) of correlation strength U. The importance 
of this kind of interference effects is known long agoi£, 
though these can be taken into account only in the case 
of weak correlations and disorder. At the same time, 
the neglect of interference effects is the key point of our 
DMFT+S approach, allowing to obtain rather simple 
and physically clear interpolation scheme, allowing to an- 
alyze the limits of strong correlations and disorder. 

Another drastic simplification is our assumption of non 
magnetic (paramagnetic) nature of the ground state of 
Anderson - Hubbard model. The importance of mag- 
netic (spin) effects in strongly correlated and disordered 
systems is obvious, as well as the importance of competi- 
tion between different kinds of magnetic ground states^. 

Despite these shortcomings, our results seems rather 
attractive and reliable, e.g. with respect to strong disor- 
der effects on Mott - Hubbard transition and the general 
form of the phase diagram at T = 0. Our predictions for 
the general behavior of dynamic (optical) conductivity 
and disorder induced Mott insulator to effective "metal" 
transition can be directly compared with existing and 
future experiments. 
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